o 
o 

(N 
Or 

c/3 



> 

On 
O 

o> 

O 

o 

S3 

Pl,: 

6 
h 

CO 



X 



The accretion history of the Universe with the SKA 

Matt J. Jarvis* and Steve Rawlings^ 

Astrophysics, Department of Physics, Keble Road, Oxford, 0X1 3RH, UK 

In this paper we investigate how the Square Kilometre Array (SKA) can aid in determining the evolutionary 
history of active galactic nuclei (AGN) from redshifts z = — > 6. Given the vast collecting area of the SKA, it will 
be sensitive to both 'radio-loud' AGN and the much more abundant 'radio-quiet' AGN, namely the radio-quiet 
quasars and their 'Type-IF counterparts, out to the highest redshifts. Not only will the SKA detect these sources 
but it will also often be able to measure their redshifts via the Hydrogen 21 cm line in emission and/or absorption. 
We construct a complete radio luminosity function (RLF) for AGN, combining the most recent determinations 
for powerf ul radio sources wi th an estimate of the RLF for radio-quiet objects using the hard X-ray luminosity 
function of Ue da et alJ 1)2003), including both Type-I and Type-II AGN. We use this complete RLF to determine 
the optimal design of the SKA for investigating the accretion history of the Universe for which it is likely to be a 
uniquely powerful instrument. 



1. INTRODUCTION 

Radio astronomy has played an important role 
in the hunt for high-redshift active galactic nu- 
clei (AGN), with very-high-redshift objects being 
pinpointed by radio pre-selection: at low radio 
frequency finding 'steep-spectrum' radio galaxies 
i Rawlings etafl Il996t Ivan Breugel et all Il999j) 
out to redshift z ~ 5; and at high frequency 
finding similarly distant 'flat-spectru m' (or 'Giga- 
Hertz -Spectrum' GPS) quasars (e.g. iHook et all 
2002). The obvious benefits of radio pre-selection 
remain, e.g. the fact that obscuration by dust 
and neutral Hydrogen is unimportant, contrast- 
ing dramatically with the case for optical and X- 
ray surveys. 

The reason why radio surveys led the way for so 
long is because a radio-loud source is essentially 
detectable out to the very highest redshifts with 
only a moderate exposure time, e.g. the highest 
redshift radio galaxy discovered to date was se- 
lected fro m the Westerbork North ern Sky Survey 
(WENSS 3 : lRengelink et alllll997fn and has a flux 
density of 5*365mhz ~ 0.6 Jy. Another advantage 
of radio surveys is that, particularly at the lower 
frequencies, the sky coverage per pointing can be 
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huge, e.g. for the new 74 MHz extension at the 
VLA the field-of-view (FOV) is ~ 40 deg 2 , much 
larger than any FOV on an optical or X-ray tele- 
scope. 

However, genuinely radio-loud sources are 
rare, probably ~ 10 — 100 times less abun- 
dant than their radi o-quiet counterparts 
dGoldschmidt et ail Il999j), w hereas the opti- 
cal (|Fanet all l2001albl l2003jh and (with the 
advent of XMM-Newton and CHANDRA) the 
X-ray surveys are now taking the lead in sur- 
veying AGN in the hi gh-redshift Universe (e.g. 
iBarger et all l2003albl) . The results of these 
surveys are providing us with a much better un- 
derstanding of both the unobscured and obscured 
AGN populations, in the case of the Sloan Di 
ital Sky Survey (SDSS; ISchneider et"afl EoQ. 
and, on the traditional ground of the radio tele- 
scope, the X-ray satellites are giving a much more 
trans parent view of the m ore obscured popula- 
tion (N orman et all 12002ft . Unfortunately, there 
are still many unresolved issues, most obviously 
highlighted by the fact that about 50% of the 
hard X-ray background is still unresolved into 
discrete AGN, the consensus being that the cur- 
rent generation of X-ray telescopes, operating 
at 0.5 — > 12 keV, are missing the highly ob- 
scured sources, with Hydrogen column densities 
N H > 10 28 m" 2 ijWilman. Fabian, k Nulserl 
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1200(1 llleda et all l200fl fCandhi et all 120041) . 

The SDSS and ultra- deep X-rays fields 
are paving the way (e .g. iBrandt et all l200li 
lAlexander et all (2003b) but the quest to deter- 
mine the nature and, possibly more importantly, 
the evolution in the obscured population is li- 
able to switch back to the realm of radio astron- 
omy with the commissioning of the next gener- 
ation of radio telescopes, most notably the Low- 
Frequency Array (LOFAR 4 ) and the Square Kilo- 
metre Array (SKA 5 ). 

The LOFAR is currently set to operate at the 
very-low-frequency (< 250 MHz) end of the ra- 
dio spectrum, where the Universe has, to date, 
not been deeply explored. This will allow an un- 
precedented view of the highest redshift AGN, 
where the intrinsic spectral shape of a radio- 
loud AGN means that the most sensitivity is 
gained at the lowest frequencies. Coupled with 
the fact that the 21 cm Hydrogen (HI) line is 
redshifted to below 200 MHz at z > 6, where 
the epoch of reionization is coming to an end 
I Becker et aTll200lh . this means that the LOFAR 
will be able to probe the neutral Hydrogen within 
the epoch of reionization, as opposed to the op- 
tical waveband where the Gunn-Peterson trough 
due to n eutral Hydrogen extingu ishes Lyman-a 
photons llOunn ,V ['etersonHl'Ki.-J). However the 
LOFAR , although important for many studies 
of the AGN population, will not be able to se- 
cure redshifts for the vast majority of the ra- 
dio sources it finds, and as such optical and 
near-infrared spectroscopy will always be needed. 
Moreover, due to the large source density at 
the faintest flux-density limits, selection crite- 
ria based on the radio spectral properties will 
still be required in order to disentangle the high- 
est redshift sources from the more prominent 
low-re dshift populatio n (e.g. Rottgering et all 
1991 ICha,mhers et all UmA fenttgering et all 
1997; Blund ell et allll998tlja7vis et all l2001afbl 



De Breuck et all l200l|) 



A more promising telescope in terms of con- 
straining the evolution and properties of the AGN 
population is the SKA. With the proposed sensi- 




Figure 1. The flux density of a 'radio-quiet' AGN 
with Li.4GHz = 10 22 W Hz" 1 sr- 1 and a = 0.8, 
as a function of redshift and observing frequency. 
The solid line is the flux-density measured at the 
observed frequency of 1400 MHz; the dashed line 
is the flux density measured at 700 MHz and the 
dotted line is the flux density measured at the 
frequency of redshifted HI. 



tivities, the SKA will be able to detect all the 
radio-loud and radio-quiet 6 AGN up to and in- 
cluding the epoch of reionization at z > 6 in 
a ~ 30- minute integration (Fig. The capa- 
bilities of the SKA for probing AGN within the 
Epoch of Reionization are considered in Carilli 
et al. (this volume). The SKA also has the ca- 
pability to measure the redshifts of many of the 
sources it detects via the HI line in emission or 
absorption. A three-dimensional picture of the 
whole of the AGN population will become appar- 
ent immediately. 

In this paper we investigate what the SKA will 
be able to contribute to studies of the whole of 
the AGN population, both 'radio loud' and 'ra- 
dio quiet' (see also Jackson, this volume; Falcke, 
Kording & Nagar, this volume). In Section |3 we 



4 http:/ /www. lofar.org 
J http:/ /www. skatelescope.org 



6 A characteristic value of Li.4GHz = 10 22 W Hz -1 sr — 1 is 
chosen for 'radio-quiet' AGN above which (roughly speak- 
ing) weak-AGN-dominated radio emitters dominate, and 
below which starburst-dominated radio emitters domi- 
nate. 
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describe the form of the radio luminosity func- 
tion (RLF) for the powerful radio-loud AGN. In 
Section |3| we show how to translate the hard-X- 
ray luminosity function to that of an RLF for the 
radio-quiet AGN population. Adding these to- 
gether yields a complete radio luminosity function 
for AGN, described in Section 0] In section 
we estimate what the SKA will see in terms of 
AGN and in Section [fj] we investigate the opti- 
mum design for the SKA in terms of constraining 
the evolution of the entire AGN population. In 
Section we outline some of the science that the 
SKA is likely to deliver. 

We adopt a ACDM cosmology with the follow- 
ing parameters: h = i/o/(100 km s -1 Mpc -1 ) = 
0.7; Q m = 0.3; fl A = 0.7, throughout. We define 
radio spectral index as S v oc where S v is the 
radio flux density at frequency v. 

2. The radio-loud luminosity function 

The radio-loud population can be split into two 
types of source, with a division (roughly speak- 
ing) at a critical radio l umino sity in radio struc- 
ture llFanaroff fc Rilevl Il974l\. narrow-emission- 
line stre ngth Jrline fc Longairl Il979t) and quasar 
fraction (Wi llott et all l2000|) . Below a radio lu- 
minosity Li.4GHz ~ 10 25 W Hz -1 sr _1 , objects 
are typically FRI radio galaxies with weak or ab- 
sent emission lines. Above the break in the RLF 
(at Li.4GHz ~ 10 26 W Hz -1 sr _1 ), objects are 
a roughly equal mix of FRII radio galaxies and 
FRII radio quasars, typically with strong emis- 
sion lines, indicating buried quasar nuclei even in 
the cases where no direct quasar emission is ev- 
ident. The transition region between these two 
sub-populations lies just below the break in the 
RLF. The vast majority of work in constrain- 
ing the evolution in the luminosity function of 
the powerful radio sources has concentrated on 
the 'FRII sub-population', due mainly to the fact 
that the low radio luminosities of FRI sources 
makes them hard to detect out to the highest red- 
shifts. 

2.1. The high-luminosity RLF: FRIIs 

The form of the RLF for the FRII sub- 
population has been constrained to a certain de- 



gree by the work of iDunlop fc Peacqckl l)l990h . 
and more recently IWillott et alJ (|200lf ). Here we 
use one of the parameterized models of Willott et 
al. (model C) which describes the high-luminosity 
population by a 'reversed Schechter function' 
with an evolutionary term htted with two one- 
sided Gaussian functions, forced to match at a 
peak redshift of z — 1.9. The form of this 
RLF means that the increase in source den- 
sity towards this redshift is much steeper than 
the hig her redshift decline. However, the data 
used by IWillott et all l|200ll) to model the RLF 
has a dearth of sources at high redshift due 
to a high flux-density limit for the large area 
samp les (e.g. 3CRR; lLaing. Rilev. fc Longairl 
Il983|) and a lack of sky area for the 
faint flux-density samples (6CE and 7CRS; 



taint flux-density sample s (bUK and /l^s; 
Fa,les et a,lllT997tlRawlm g s. Kales, fc LacvlEoOll 

IT 



Willott et all bond I2003D . Thus, the high- 
redshift evolution is not yet well constrained. 

The form of the decline in the comoving space 
density at high redshift redshift has come under 
close scrutiny from a number of groups, some 
of which suggest qu ite a steep decline at z > 3 
^Shaver et all Il996j) . and so me of which suggest 
a more moderate decline (e.g.lDunlop fc Peacockl 

e . However Ijarvis fc Rawlingj <)2000|) and 
3 et all I 2001c^ highlighted how uncertain 
the form of this decline is for both flat-spectrum 
radio-loud quasars (which, crudely speaking, 
should reflect the overall evolution as it picks 
out favourably oriented sources from the under- 
lying population) and also low-frequency-selected 
radio galaxies, due predominantly to difficulties 
in modelling the k— corrections, and the lack 
of available volume (and hence small numbers 
of sources) at high redshift for a given flux- 
density-limited sample. This uncertainty is also 
mirrored by the X-ray luminosity function (e.g . 
[ Mivaii. Hasinger. fc Schmidtl 120001 llleda et all 
2003, see also Section |3J|. Our adopted model lies 
between the extremes of a sharp cut off and an 
RLF which stays constant with redshift at z > 2. 
Such an evolution is consistent with the predic- 
tions of simple Press-Schechter-b ased theoretical 
models ifoawlings fc Jarvisi Eool . 
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2.2. The moderate-luminosity RLF: FRIs 

Unfortunately, the RLF for the lower- 
luminosity radio-loud AGN, predominantly (but 
not exclusively) comprised of FRI radio sources, 
hereafter the 'FRI sub-population', is poorly 
constrained. Several groups have reached 
different conclusions 
tionary behavi our (e.. 



concerning their 


evolu- 


g. Jackson & Wall, 


1999; 


IClewlev & Jarvisl 12004). 



(2004), who used the largest dataset of ~ 
1000 sources, with a radio flux-density limit of 
S325MHZ = 40 mjy, with photometric redshifts de- 
termined by five-colour photometry of the SDSS, 
to show that radio sources with 10 23 < (£325MHZ 
/ W Hz -1 sr^ 1 ) < 10 25 are consistent with 
having a constant comoving density, at least in 
the redshift range z — — > 0.5. We also 
note that the res ults of their analysis are c onsis- 
tent with that of Waddin gton et alJ (|200 If) . who 
constrained the form of the evolution using a 
much fainter flux-density limited survey (LBDS; 
IWaddington et all 120001 but with a higher frac- 
tion of redshift incompleteness. We allow this 
sub-population to have a comoving space den- 
sity which evolves very mildly with redshift as 



$ cx (1 



\1.2 



, normalised to the comoving space 



density of this sub-population from the model of 
Willott et al. at z = 0.5, i.e. where redshift data 
are able to directly constrain the space density. 

3. The low-luminosity RLF: radio-quiets 

As stated in Section the SKA will have the 
sensitivity to trace radio-quiet quasars and radio- 
quiet type-II AGN out to the highest redshifts. 
Therefore, we have attempted to parameterize 
the form of the radio-quiet RLF and its evolu- 
tion w ith redshift. The work of iKukula et al.l 
l|l998f ) suggested that the radio emission from 
radio-quiet quasars is typically compact and is 
probably directly associated with the central en- 
gine. They also showed that radio-quiet quasars 
have spectral indices consistent with their more 
radio-loud counterparts, i.e. a ~ 0.7. 

We also know that the X-ray luminosity of an 
AGN is a good tracer of the accretion rate onto 
the SMBH. Thus, there should be a link between 



the X-ray luminosity and the radio luminosity if 
both are regulated by the acc retion process. In- 
deed, [Brinkman^^^O l|2000l) show that this link 
does exist. They compiled a sample of X-ray 
selected AGN from the ROSAT All-Sky Survey 
l|Voges et all l"l999) and cross-correl ated this sam- 
ple w ith the 1.4 GHz FIRST survey (|White et all 
Il997l) . They find a relatively tight linear corre- 
lation between the logarithm of radio luminosity 
at 1.4 GHz and the X-ray luminosity in the soft 
X-ray band at 0.5-2 keV band. The form of this 
correlation is 



log(ix) = -4.57+1.012 log(Li 



(1) 



i.e. roughly a proportionality. With this in mind, 
it is possible to use a complete description of the 
radio-quiet quasar population in the X-ray band 
to make a crude estimate of the RLF of radio- 
quiet objects at 1.4 GHz. 

With the advent of XMM-Newton and Chan- 
dra, and the deep fields surveyed with these satel- 
lites, a hard X-ray lu minosity function (XLF) 
has been determine d bvllJeda et alJ ll2003lL build- 
ing on the work of iMivaii. Hasinger. fc Schmidtl 
(2000) in the softer X-ray bands. They find 
that the best parameterisation of the XLF is 
one of luminosity-dependent density evolution, 
i.e. one in which the cosmic evolution in the 
space density of the sources depends on the in- 
trinsic luminosity of the source. They conclude 
that the comoving space density of the more X- 
ray-luminous AGN peaked at an earlier epoch 
than those at lower X-ray luminosities. The ev- 
idence for this ty pe of evolution is re inforced by 
other studies re.g. lSteffen et all 120031) . They also 
find that the higher luminosity sources have less 
intrinsic absorption, consistent with a 'receding- 
torus-like' model (e.s. lLawrenceL fl99ll) which has 
been suggested as a necessary physical ingredi- 
ent to expla in the propert i es of radio-selected 
samples (e.g.lSimpsoilll99StlWillott et all l2000: 
iGrimes. Rawlings. fc Willottl l2004|) ~ 

With this parameterization, we predict how 
many X-ray sources there will be of a given lu- 
minosity in the (rest frame) 2 — > 10 keV en- 
ergy range, at a given redshift. We use the rest- 
frame ratio of the average number of photons in 
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this band to that of the monochromatic lumi- 
nosity at 1 keV for an unobscured source, not- 
ing that the quasars were predominantly unob - 
scured in the study of iBrinkmann et alJ l)2000|) . 
We then convert this soft-X-ray luminosity to the 
1.4 GHz radio luminosity using Equation^ In 
addition, we need to take some account of a popu- 
lation of highly-obscured (Compton thick) AGN, 
which, whilst essentially absent from the exist- 
ing 'hard' (2-10 keV) surveys , would be revealed 
by st ill harder X-ray surveys IjWilman fe Fabiari 



1999) with the next generation of X-ray satellites 
(e.g. XEUS 7 and Constellation-X 8 ) or radio sur- 
veys. To do this we multiply the comoving space 
density derived from the XLF by a conservative 
factor of 1.5 to account for the X-ray obscured 
sources, which brings the integrated luminosity 
closer to that expected from the hard- X-ray back- 
ground. 

4. The complete RLF and constraints from 
the source counts 

In Figure we show the complete RLF, in- 
corporating the radio-loud sources (consisting of 
the two FRI and FR II sub-populations following 
IWillott et al.l I^OPl]) ) and the radio-quiet quasars. 
This is the RLF we use to predict what the SKA 
will see in terms of AGN for various SKA designs. 
With this parameterisation of the RLF we can 
make a consistency check with the radio source 
counts at 1.4 GHz. In Figurc[2]we show the radio 
source counts at 1.4 GHz down to a flux-density 
limit of Si.4GHz ~ 1 /uJy, along with the source 
counts derived from our RLF. 

One can see immediately that the RLF 
traces the source counts reasonably well at the 
high-flux-density regime, as it should do, as 
this is where direct constraints from data ex- 
ist. The small deviations are artefacts of the 
'two-sub-population' method which would be 
removed by using smoothe r functional forms 
(such as the 'free- form' fits of lDunlop fe Peacock! 
lll990l) or the generalized luminosit y functions of 
lOrimes. Rawlrngs. fc Willotil (|£o04) . At lower ra- 
dio flux densities the moderate-luminosity 'FRF 



7 http:/ /www.rssd.esa.int/XEUS 
8 http: / /constellation. gsfc.nasa.gov 



AGN begin to dominate the source counts, but at 
lower flux densities still (Si.4gh z ~5x 1CT 4 Jy) 
the radio-quiet population can start to contribute 
significantly to the source counts. The upturn 
in the source counts around this flux density is 
often attributed solely to a low-redshift, star- 
forming population. However as can be seen 
from Fig. the radio-quiet population, compris- 
ing radio-quiet quasars and their 'Type-II' coun- 
terparts (including Compton-thick sources) may 
contribute significantly to this upturn (see also 
Baugh et al., this volume). Of course, starburst- 
driven radio emission may make an important 
contribution to the radio flux densities of at least 
some radio quiet quasars ( e.g. lSopp fc Alexanderl 
ll992HProuten et alJ.E004^ . and 'starburst' galax- 
ies seem commonly to harbour X -ray (AGN) nu- 
clei ('e.g. I Alexander et all l20f)3a|) . so the distinc- 
tion between these two populations is not par- 
ticularly clear. The starburst population is con- 
sidered elsewhere (van der Hulst et al., this vol- 
ume) , and it is clear from Fig. |3| that as we ap- 
proach the depths of even shallow SKA surveys 
(surveys limited at ~ /iJy levels), a small frac- 
tion of the objects are likely to have their radio 
emission powered by any sort of AGN activity. 

Our description of the radio luminosity func- 
tion, although very crude, allows us to at least 
estimate the contributions from the various pop- 
ulations of radio AGN. As such we can distin- 
guish between those with massive black holes but 
low accretion rates (FRIs) and those with mas- 
sive black holes and high accretion rates (the 
radio-qu iet quasars). However, it i s worth not- 
ing that Blu ndell fc Rawlinga (|2001f ) have shown 
that both classes can produce similar radio lumi- 
nosities and radio structures. Ideally, a complete 
model for the evolution of all AGN in radio wave- 
bands would allow for more physically driven pro- 
cesses, with a more continuous level of accretion 
between the various classes (or equivalently less 
of a dichotomy between 'radio-loud' and 'radio- 
quiet' sources). 

Also, we have not attempted to add prescrip- 
tions for favourably-oriented, and hence Doppler- 
boosted (flat-spectrum) sources and young (GPS) 
sources as we will detect every AGN with the 
SKA and these sources do not make up the bulk 
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Figure 2. The 1.4 GHz radio luminosity function for all AGN at four different redshifts. The contributions 
from the radio-quiet quasar population can be seen to be significant between 10 22 < (X1.4/W Hz -1 sr^ 1 ) < 
10 24 , whereas the bulk of the radio luminosity density is provided by the FRI-type sources at L1.4 ~ 
10 25 W Hz -1 sr -1 . Finally, at L1.4 > 10 26 W Hz -1 sr -1 all of the luminosity density is provided by the 
powerful FRII-type sources. The thin solid line is the RLF for each sub-population at z = 0, the dotted 
line is at z = 0.5, the dashed line is the RLF at z = 1 and the dot-dashed line is the RLF at z = 2. To 
clarify the overall form, the thick solid line represents the composite RLF, comprised of the FRII, FRI 
and radio-quiet populations described in Section |2 at z = 2. 



of the population when probing such low flux- 
density limits. 

Again, more realistic models could allow for 
factors such as the time evolution of radio lu- 
minosity (although we caution that such models 
bring in lots of uncertain physical details such as 
evolutions in environment, and the behaviour of 
jet power as a function of time). 

However, regardless of all such details, with a 
crude composite RLF in place we are able to make 
a crude prediction of the number of each type of 
AGN the SKA will be able to detect as a function 
of flux density and area. 



5. What will the SKA see? 

5.1. The number of AGN detected by the 
SKA in 1 year 

Using the composite RLF discussed in Section|21 
we are able to estimate the number of sources 
that the SKA will be able to detect as a func- 
tion of both redshift and flux-density limit. In 
Fig. 0] we plot the predicted dN/dz distribution 
for each sub-population in an all-hemisphere SKA 
survey which would reach a limiting (4cr) flux- 
density limit of 52oomhz = 100 nJy in one year of 
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Figure 3. The 1.4 GHz differential source counts (re l ative to counts in a Euclidean Universe) of a 
combination of surveys fsee lSevmour. McHardv fc Gunrl 120041 for a description of the data points). The 
lines represent the model source counts from our three components of the radio luminosity function. The 
dotted line represents the powerful FRII sources, the dot-dashed line is for the moderate luminosity, 
typically FRI sources, and the dashed line represents the source counts for radio-quiet quasars. The solid 
line is the total source counts from these models. 



dedicated observations with the SKA 9 . 

One can see that the SKA will be able to de- 
tect every bona fide AGN in the Universe up to 
z ~ 6. It is worth noting that there is a large over- 
lap between what we have called FRI-type AGN 
and the radio-quiet quasars. This is because it is 

9 This assumes a usable field of view at 200 MHz of 
« 50 deg 2 (or FOV = 1 deg 2 at 1.4 GHz in the diffraction- 
limited case), so a survey of ~ 20, 000 deg 2 requires ~ 400 
pointings. A one-year survey (at 90 per cent observing effi- 
ciency), means ?a20 hours of exposure per pointing, yield- 
ing a survey rms of 25 njy. In practice, a 'tiling' method 
will be used to ensure complete sky coverage across the 
entire 0.2 — » 1.4 GHz range (see Rawlings et al., this vol- 
ume). For an SKA design with a ~10-timcs-higher FOV, 
a 10-it catalogue to the same limiting flux density could 
be compiled from a survey taking eight months. 



largely a question of seman tics as to what kind of 
sourc e belongs where (e.g. iBlundell fc Rawlingst 
120011) . 

It is also worth noting that at high redshift 
(z > 3) all three of the sub-populations are 
largely unconstrained. Therefore, we force both 
the FRI-type sources and the radio-quiet sources 
to decline according to the decrease in space 
density expecte d under a Press-Schechte r based 
model (see e.g. iRawlines fc Jarvisl l2004j) above 
a redshift of z ~ 3. For the more powerful, 
FRII-type sources we just use the shallow de- 
cli ne in the space den sity described by model C 
in lWillott et~ai1 l|200l[) . 

Due to these large extrapolations the ratio be- 
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Figure 4. The number of the different sub-classes 
of sources discussed in Section 0] for a one year 
SKA survey over one hemisphere (~ 20, 000 deg 2 ) 
complete to S200MHZ ~ 100 n Jy- The dotted line 
represents the powerful FRII-type sources, the 
dashed line shows the evolution in the FRI-type 
sources and the solid line denotes the radio-quiet 
quasa rs calculated using the XLF of lUeda et aD 
(2003) and multiplied by a factor of 1.5 (see Sec- 
tion EJ. 



tween radio-loud sources and radio-quiet sources 
stays roughly similar to what we observe at z < 2 
with today's telescopes. However, the SKA will 
actually measure this ratio to great accuracy, and 
as a consequence may yield strong constraints on 
AGN driven feedback ev ents (see Section and 
iRawlings fe Jarvisl l)2004[) ') and the duty cycle of 
all types of AGN. 

5.2. Supermassive black holes and their ac- 
cretion rate 

In this section we estimate what the SKA will 
be able to see in terms of the accretion rate of 
black holes with masses in the range 10 6 M < 
Afbh < 10 9 Mq. To estimate the bolometric lu- 
minosity, and thus calculate the fraction of the 
Eddington luminosity which relates to a given ra- 
dio luminosity for both the powerful radio-loud 
sources and the radio-quiet sources we use two 



separate methods. 

For the radio-loud sources, i.e. those sources 
above the brea k in the RLF (see Se ctionEt we use 
the relation of lWillott et alJ l)2000|) which relates 
the low-frequency radio luminosity to the total 
ionizing power of the AGN. This has the form, 



]og(£ Bo j/W) « 0.831og(L 15 i/W Hz^sr"^^.^ (2) 



where LboI is the bolometric luminosity and L151 
is the monochromatic luminosity at 151 MHz. We 
convert the 151 MHz luminosity to 200 MHz using 
a spectral index of a = 0.8. 

For the case of the radio-quiet quasars (and 
type-II obscured quasars) we use the ratio of bolo- 
metric luminosity to monochromatic lum inosity 
L 250 o (at 2500 A) from lElvis et alJ l)l994[) - along 
with the newly deriv ed optical — X-ray spectral 
index «ox given in lElvis. Risaliti fc Zamoranil 
(2002J). We then use the correlation between 
X-ray luminosity and radio luminosity given in 
Section |31 to estimate t he radio lum i nosity at 
200 MHz (using a = 0.7: iKukula et allll99s() . 

Fig-GDshows what the SKA, with a flux-density 
limit of 52oomhz > 100 nJy, will see in terms of 
radio-loud quasars (left) and radio-quiet quasars 
(right) as a fraction of the Eddington luminosity 
versus redshift for AGN with black-hole masses 
10 6 M Q < Mbh < 10 9 M Q . Given that power- 
ful radio sources seem to require a supermassive 
black hole with mass Z 1 8 Mk v, (iDuiil^c^g^lJ, 
| 2003[ I.Tarvis fc McLureL 12003: iMcLure et all 
200l IMcLure fc Jarvisl Eool then it is obvi- 
ous from Fig. [5] that the SKA will easily de- 
tect all such sources. However, more interest- 
ingly the SKA will be able to detect all the radio- 
quiet quasars with Mbh > 10 6 M Q accreting at 
< 10 per cent of its Eddington limit up to z ~ 6. 
In the case of the powerful quasars being found 
at high redshift in t he SPSS, with black-hole 
mass es > 10 9 M l|Willott. McLure fc Jarvisl 
2003), these would be detected up to z ~ 6 if their 
accretion luminosities were ~ 10 -4 L^,dd, i-e. in 
the realm of what we now consider as dormant 
black holes. 
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Figure 5. The detectable accretion luminosity L/L^dd versus redshift for sources containing black holes 
of various mass (10 6 — > 10 9 M Q ), at a flux-density limit of S200MHZ = 100 nJy for radio-loud, typically 
FRII, sources (left) and for radio-quiet sources (right). 



6. The optimal SKA design 

In this section we investigate the most efficient 
design of the SKA for determining the large scale 
properties and evolution of the AGN. We only 
consider those AGN which have an X-ray lumi- 
nosity Lx > 10 35 W, i.e. those in the regime 
where accretion onto a supermassive black hole 
is the only way of producing such emission and 
avoiding luminosities where supernovae and X- 
ray binaries may begin to dominate the popula- 
tion. This means that the SKA will be able to 
probe the 'whole AGN population' up to z ~ 6 
by just surveying down to micro-Jansky levels. 

One of the crucial points in a survey to map 
the AGN population at < z < 6 with the SKA 
is that IT MUST be able to observe across the 
frequency range 200 < v < 1400 MHz. This is 
not for the traditional reason of previous low- 
frequency-based redshift surveys such as 6CE, 
6C* and 7CRS. In these samples the low sur- 
vey frequency selected sources in an unbiased 
manner with respect to source orientation, i.e. 
sources are selected on their optically-thin lobe 
emission. Low-frequency selection is also used to 
target the hi g hest redshift radio galaxie s (see e.g. 
ICohen et all 12004 Ijarvis et al.ri2004j) . whereas 



traditional high-frequency selection would miss 
them due to the shape of radio galaxy spectra. 
The sensitivity of the SKA will mean that this 
type of 'low-frequency selection' is unimportant, 
as all of the radio sources will easily be detected 
at 1.4 GHz regardless of the redshift. 

However, the crucial advantage that the SKA 
will have over future X-ray missions in probing 
the AGN population will be its ability to mea- 
sure re dshifts of all the gas-rich g alaxies out to 
z ~ 2 l|Abdalla fc Rawlingd. |2004|) via HI emis- 
sion. This means that the redshifts of the ma- 
jority of the source counts at these flux-density 
levels will already be known. This paves the way 
for both directly constraining the evolution in the 
high-redshift population along with their three- 
dimensional clustering, and to investigate the ac- 
cretion rate of all powerful AGN from z — — * 2. 

Moreover, powerful AGN, invariably reside in 
massive galax i es (i.e . > L*: Ijarvis et all l2001ri 
IWillott et all l2003t iMcLure et all l2004|) and 
high-redshift- AGN activity, triggered by mergers, 
are likely to contain higher gas fractions than at 
lower redshifts. Evidence f or this has been seen 
in the past few years (e.g. IPanadopoulos et all 
l200ll iGreve. Ivison k Papadopoulosl l2004|) 

Therefore, the HI emission signature will prob- 
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ably be detected to higher redshifts than for 
the less massive gal a xies w hich are detailed in 
lAbdalla k Rawlingj l|2004|) . There is also the 



possibility of utilizing the fact that many AGN, 
both optically selected and radio selected, exhibit 
some level of HI absorption in their Lyman-q 



profiles, (e.g. 


Ivan Oiik et ; 


,11119971 I.Tarvis et all 


2003; 


Wilman et all 


2004 


). This gas could also 



be detected against the radio continuum of these 
distant sources via the HI transition. However, 
this obviously needs the SKA to operate at lower 
and lower frequencies depending on the maximum 
redshift one wishes to probe. 



In order to estimate whether the SKA will 
be able to determine redshifts via HI absorption 
against radio-loud objects up to z ~ 6 (see also 
Morganti et al., this volume) we assume a min- 
imum HI column density of Ahi = 10 23 m~ 2 
(or AThi = 10 19 cm~ 2 ) and a line width of ~ 
100 km s" 1 which is similar to the absorption sys- 
tems observed in Lyman-q around high-red shift 
radio galaxies, ijjarvis et alll2003l:IWilman et all 
l2004) . 

Usin g equation 1 from IChengalur fc Kanekarl 
i|2000f) w ith a spin temperatu re of T s = 100 K 
(see e.g. iBeswick et all 12004) and unity cover- 
ing factor we find a target velocity-integrated 
optical depth rdV ~ 0.06 km s _1 . Taking a 
source at z ~ 6 with a luminosity around the 
break in the RLF for powerful radio sources, i.e. 
£200 ~ 10 26 - 5 W Hz -1 sr _1 gives a flux density 
at the redshifted frequency of HI of S^oomhz ~ 
15 mJy In the proposed one- year all-sky survey, 
the expected dip due to HI absorption (at a level 
~ 6 x 10~ 4 of the continuum) would be detected 
at w 10(7 in an optimally-smoothed spectrum. 10 

Obviously, ignoring any cosmic evolution in the 
HI properties, the situation 11 becomes much bet- 
ter at lower redshift where higher sensitivities are 
reached for a source of the same luminosity as 
the one considered at z ~ 6. Therefore, we en- 
visage that the SKA will be able to detect and 
measure the redshift of almost every 'FRII-like' 
AGN in the Universe up to at least z ~ 6, if 

10 Assuming FOV = 1 deg 2 . If FOV was (say) 10-times 
higher, these features would be detected at ss 30<r and 
weaker (lower TVjj and higher Tg) features could be de- 
tected. 



these sources lie within a shell of HI gas with col- 
umn densities ~ 10 19 cm -2 , as has been observed 
around young powerful sources at high redshift 
via resonant HI absorption within the Lyman-q 
emiss ion line ijjarvis et all l2003t IWilman et all 
l2004h . Redshifts of radio sources with lower HI 
column densities and / or higher spin temperatures 
could also be measured if their intrinsic linewidths 
were narrower, or if their continuum sources were 
brighter (see also Kanekar & Briggs, this volume). 
With HI emission detectable at low redshift, and 
the hope that HI absorption is close to ubiqui- 
tous at high redshift, an almost complete census 
of the evolution of powerful radio AGN activity 
may well be obtained using just the SKA. 

7. Conclusions 

In this paper we concentrate on z < 6 AGN 
(for a discussion of higher-redshift AGN, see Car- 
illi et al., this volume). Although detecting radio 
emission from z < 6 AGN will be a trivial exercise 
for the SKA, measuring their redshifts and prob- 
ing them in detail requires low-frequency capabil- 
ities. A SKA which can operate between 200 MHz 
— * 1400 MHz would therefore be ideal for investi- 
gating AGN activity after the epoch of reioniza- 
tion. Here we highlight some of the science that 
the SKA could deliver regarding this population. 
We should be able to achieve the following. 

• Detect all the > 10 6 M Q black holes in 
the observable Universe with accretion lu- 
minosities above ~ 0.1 ^Edd- 

• Trace accretion luminosities from 
L/Lftdd = 1 — > 10~ 4 as a function of 
redshift up to z — 6 for all black holes 
with M > 10 9 M . The SKA will also be 
able to probe the M > 10 s M dormant 
(L/I/Edd < 10 -3 ) black hole population out 
to z ~ 6. 

• Constrain the evolution in the radio lu- 
minosity function for both radio-quiet and 

To obtain full sky coverage for HI at all redshifts in the 
range < z < 6 the continuum survey will need to adopt 
the 'tiling' approach discussed in Sec. 2 of Rawlings et al. 
(this volume). 
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radio-loud AGN, up to z ~ 2, without the 
need for optical follow-up time. This is pos- 
sible because the SKA will be able to de- 
tect HI emission from all massive galaxies 
up to this redshift for a survey in which the 
effective exposure time is ~ 4 — 20 hours 
(see Rawlings et al., this volume). Con- 
straining the luminosity function at higher 
rcdshifts, i.e. within the 'quasar epoch' at 
z ~ 2 — > 6, requires (if additional opti- 
cal or near-infrared observations are to be 
avoided) HI absorption experiments with an 
SKA that operates down to 200 MHz (see 
also Kanekar & Briggs, this volume). 

• Constrain the three-dimensional clustering 
of these sources, with respect to the general 
galaxy population (see van der Hulst et al., 
this volume), again without the need for op- 
tical follow-up time. This may be crucial to 
our understanding of how AGN activity is 
triggered with respect to the large-scale en- 
vironment, i.e. where does the star forma- 
tion occur in a cluster and where does the 
AGN activity occur. Does it occur along 
dark-matter filaments with AGN at their 
centres? Progress will have been made on 
this issue by the time the SKA starts to 
operate, but only the SKA will trace this 
evolution in an 'all hemisphere' sample up 
to z ~ 2. 

• The huge cosmic volumes (in comparison 
to present and future optical and X-ray 
telescopes) surveyed by the SKA would 
allow investigations o f the largest struc- 
tures in the Universe. iBrand et alJ l)2003|) 
have shown that moderately powerful ra- 
dio galaxies are able to trace structures on 
~ 100-Mpc scales. On discovering such 
structures the SKA will be able to probe 
where the AGN and star-formation activ- 
ity occurs with respect to the large-scale 
structure. Counts of these structures have 
cosmological implications (Blake et al., this 
volume; Rawlings et al., this volume). 

• The sensitivity of the SKA to the HI line 
will allow extremely detailed investigations 



of the neutral gas surrounding the cen- 
tral engine, especially with VLBI capabil- 
ity (Morganti et al., this volume). This 
high column density gas (Ah > 10 26 m -2 ), 
may be responsible for the archetypal type- 
II quasars and could easily be detected with 
the SKA for any source with a moderately 
bright core, up to the highest redshifts. Not 
only will the SKA be able to detect it but 
it will also be able t o trace the dynamics 
l)Beswick et all l2004|) . This sort of study 
will never be achieved with optical tele- 
scopes because the high-column densities 
involved lead to a heavily saturated absorp- 
tion profile in Lyman-a. 

• Measure inflows/outflows of HI gas, both 
in emission and absorption up to z ~ 
2. The SKA would also have the 
power to see gas stripping by/around an 
AGN by tracing the Hl-emitting (e.g. 
iKennev. van Gorkum fc Vollmeii 12004^ . If 
powerful AGN lie at the interface of merg- 
ing s ub-clusters (e.g. ISimpson k. Rawlingst 
l2002j) . then the gas stripping around galax- 
ies due to their interaction with the intra- 
cluster medium should be easily detectable 
with the SKA. Outflows of gas, driven by 
the AGN, may strip gas from not only the 
host galaxy but also neighbouring galaxies. 
This is one area where a multi-wavelength 
approach, particularly with optical/near- 
infrared telescopes, would enable the inves- 
tigation of gas striping along with chemical 
enrichment within the intracluster medium, 
and how this is affected by the geometry of 
any sub-cluster merger. 

• An 'all-hemisphere' survey, deep enough to 
detect all the AGN in the Universe, and 
measure redshifts for all of the most power- 
ful ones, will also provide the perfect sam- 
ple for probing the magnetic fields along the 
line of sight via Faraday Rotation measure- 
ments (Gaensler, Beck & Feretti, this vol- 
ume). As the SKA will detect AGN at all 
redshifts then any evolution in the magnetic 
fields could also be investigated with such a 
survey. 
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• The SKA would greatly benefit investi- 
gations of feedback mechanisms in the 
high-redshift Universe where AGN-driven 
feedback events may be a crucial fac- 
tor in the evoluti o n of all galaxies. In 
iRawlings fc Jarvisl (|2004f) we show that 
powerful radio jets may have a profound 
influence on the evolution of all galaxies 
surrounding the AGN. This is because the 
radio source is able to inject enough en- 
ergy into its surroundings that it gravita- 
tionally unbinds ionized gas associated not 
only with its host galaxy, but more widely 
throughout the pr otoclusters in wh i ch the y 
seem to reside (jVenemans et all [2002). 
This process would yield a reservoir of 
gas within the protocluster which is grav- 
itationally unbound, and has been heated 
such that it cannot accrete back onto the 
protogalaxies. Hence, there will be a 
protocluster-wide shut down of activity, 
whether it be circumnuclear star formation 
or black-hole accretion. With the SKA we 
will be able to trace the neutral HI gas in 
such protoclusters and trace the effect that 
the powerful radio jets have on the proto- 
cluster environment. This will essentially 
be monitoring AGN feedback mechanisms 
in action, as the SKA will have both the 
sensitivity and the field-of-view to do this 
for all powerful AGN at all redshifts. 
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